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Note 


Cross-references in this handbook are indicated by the use of CAPITALS. 
In addition, references are given to appearances in the course of 
definitions/formulas/theorems. The set books are denoted by S and Y, 


and the correspondence text of Unit n is indicated by n. 


] GLOSSARY 


ABSOLUTELY CONVERGENT 
6: page 13 

l4: page 7 
AMPLIFICATION MATRIX 


15: page 7 


ANGULAR FREQUENCY 


7: page 12 


ANGULAR MOMENTUM 


16: page 13 


ANGULAR VELOCITY 


16: page 12 


oo 


The series 2 a, is absolutely convergent if 
n=1 
oo 


the series X EN converges. 
n=1 


The THREE-LEVEL SCHEME Au. = Bu. + Cu. is 
“+1 ^3 =j-1 
equivalent to a pair of two-level schemes 


which may be written in the form 


“jall ü; 
=P 
Me "à 
where 
lada Alo 
P= 
I 0 


P is called the amplification matrix. 
One-dimensional wave motion which is described 
by Asin(kx-Wt), for example, where x represents 
the displacement and t the time, has angular 
frequency W, wave length 2T/k and wave number k. 
The angular momentum about the point P of a 


particle at the point Q with velocity y and mass m 4 


mr x y, where r = PQ 
The angular velocity about the point P of a 
particle at the point Q with velocity Y is 
(xxy) /x? where r = PQ; it measures the rate of 


rotation of the particle about Q. 
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ASYMPTOTIC RATE OF 
CONVERGENCE 


ll: page 20 


BAND MATRIX 


11: page 16 


BAND WIDTH 


BESSEL'S EQUATION 
14: page 6 


W : page 179 


BESSEL FUNCTION 


Y : page 179 


Handbook 


The asymptotic rate of convergence of the 
ITERATIVE SCHEME 


xt) _ 


is -log p, where p is the SPECTRAL RADIUS of G. 


A band matrix is one in which all the nonzero 
entries are confined to the leading diagonal and 
those diagonals close to it. If ais is that 
element of a matrix A which lies in row i and 
column j, then A is a band matrix of band width 


2k+1 provided 
a.. = 0 for li-jl > k. 


1] 


See BAND MATRIX. 


Bessel's equation of order m is 


2 
ete Be (t? - m2)u = 0 
dt 
or 
d du m? 
a x) quttu=0 


A solution of BESSEL'S EQUATION (which is bounded 
at the origin) is called a Bessel function (of the 
first kind). A Bessel function of the first kind 
of order m (>0) is given by 

m+2k 


“o k 
z CD Gt) i 


eer = weg Klm) 


BODY FORCE 


l : page 10 


BOUNDARY CONDITIONS 


l: page 9 


BOUNDARY VALUE PROBLEM 
CANONICAL FORM 


CAUCHY PROBLEM 


2: page 18 


CENTRAL-DIFFERENCE 
OPERATOR 


3 page 9 


CHARACTERISTICS 


Glossary 7 


A body force is a force which acts on each 
element of a body from outside the body (e.g., 


gravitational force). 


Boundary conditions are conditions at a spatial 


boundary for all times. For example, the 


equation 
3? 9? 
a "gu x € (0,8), t > 0, 


may be given with the boundary conditions 


Vv 


u(a,t) = f(t) t20, 


3*(g, c) = g(t) t 


iv 


0. 
See GREEN'S FUNCTION. 
See STANDARD FORM. 


A Cauchy problem (or initial-boundary value 
problem) is a problem where the solution 

to a partial differential equation is required 
subject to the solution and its normal 
derivative both being specified on the initial 


line. 


The central-difference operator 9n is defined by 


au : xpeyu(xtth) - u(x-]h) xe R. 


See CLASSIFICATION OF OPERATORS. 


CLASSIFICATION OF 


OPERATORS 
2: page 19 
W: page 43 


COEFFICIENT OF 


VISCOSITY 


COMPATIBLE 


COMPLETE 


COMPLEX FOURIER 
SERIES 

6: page 18 
CONDITIONALLY 
STABLE 


8: page 17 


CONSISTENT 


8: page 9 


PDE Handbook 


The operator 


3?u MER du du duy; 
L:u > adem *BO Ot KENT +F k tun pra is 


(a) hyperbolic at (x,t) if [B(x,t)]?- 4A(x,t)C(x,t) 


v 


0, 


(b) parabolic at (x,t) if [B(x,t)]?- 4A(x,t)C(x,t) = 0, 


^ 


(c) elliptic at (x,t) if [B(x,t)]?- 4A(x,t)C(x,t) < 0. 
For hyperbolic equations, the two families of curves 


which satisfy the equation 


are called characteristics. 


See VISCOSITY. 


See CONSISTENT. 
See CONVERGENCE IN THE MEAN. 


The FOURIER SERIES expansion of a function (with 
domain R or C) in terms of the ORTHOCONAL BASIS 
le: ne Z, x e [-mn,n]) 


is a complex Fourier series. 


A FINITE-DIFFERENCE SCHEME is conditionally stable 


if it is STABLE under certain conditions (e.g., for 
particular values of the MESH RATIO) and unstable 
under others. 

A FINITE-DIFFERENCE REPLACEMENT of a differential 
equation is consistent (or compatible) if its 

LOCAL TRUNCATION ERROR approaches zero as the MESH 


spacings tend to zero. (See also LAX'S THEOREM.) 


CONSISTENT ORDERING 


11: page 28 


CONTINUOUS WITH 
RESPECT TO DATA 
W: page 6 
CONTINUOUSLY 


DIFFERENTIABLE 


Glossary 


In the successive computation of a set of 

values {uje sug} using the GAUSS-SEIDEL or 

SOR iterative scheme, a scanning order of the 
points 1,2,...,N is a consistent ordering if it 
yields the same set of equations (at each 
iteration) as would be obtained from some matrix 


which can be partitioned as 


Y E 
E D E, 
Baa Pmt Fm 
E D 
m-1 m 


where the D; are diagonal matrices and m 2 2. 


A problem is continuous with respect to its data 


if solutions corresponding to data which differ 


by small amounts also differ by small amounts. 


A function f is continuously differentiable on 


the interval I if its derived function f' is 


continuous on I. 


CONVERGENCE 
IN THE MEAN 


W: page 71 


CONVERGENCE OF 4 
FINITE-DIFFERENCE 
REPLACEMENT 


8: page 8 


PDE Handbook 


The infinite sequence of functions SEERE 


converges in the mean to f (on the interval |a,b]l) 


with respect to the positive weight function p if 
lim [oce-sy)?e = 0% 
Nao 
Let (6) be an infinite sequence of functions 
ORTHOGONAL on [a,b] with respect to the weight 
function p. For any function f let Sy = Z i= 1 
be the sum of the first N terms of the FOURIER 
SERIES of f. If the limit above is zero for every 


function f for which fiso is finite, the set lo? 


is said to be complete. 


A FINITE-DIFFERENCE REPLACEMENT of a differential 
equation is convergent if, as the MESH spacings 
tend to zero, the finite-difference solution tends 
to the true solution of the differential equation 
either at a fixed point, or for all points along 


the furthest time level under consideration. 


A subset of R? (or R?) is convex if no straight 
line intersects its boundary in more than two 
points, except possibly along a straight-line 


(or plane) section of the boundary. 
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CURVE A curve C in R? is given by 


iN 


page 18 (Gy): x = $(t),y = 00) t <t<t} 


i 


page 49 
where $ and O are continuous. If $ and 6 are 


CONTINUOUSLY DIFFERENTIABLE on [ToT] with 

0'? + $'? > 0 we say that C is a continuously 
differentiable curve. If these conditions hold 
except for a finite number of values of T, Cis a 


piecewise continuously differentiable curve. The 


curve C is closed if 


$0) - $(T)) and 8) = 6g). 


D'ALEMBERT'S SOLUTION  D'Alembert's solution is the general solution of 


W: pages 9,13 the one-dimensional wave equation 


Bull du 


9x? c? ar? 
in the form u(x,t) = p(xtct) + q(x-ct), where 


p and q are arbitrary functions. 


DELTA FUNCTION The delta function in n dimensions is defined by 
10: page 15 &(r-p) = 0 r*p 
f(p peD 
[Esper x: 
f 0 péD 


for all functions f and all domains D c R". 


DIFFUSION EQUATION The diffusion equation is 
3: page 31 Ou _ KV2u = 0 
9t 


DIRECTIONAL DERIVATIVE The directional derivative of the SCALAR FIELD 
3: page 7 $ in the direction of the unit vector e is given 


by the inner product e*grad. 
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DIRICHLET PROBLEM A boundary value problem in which a solution, u, 
3: page 9 of a differential equation in a DOMAIN is 
9: page 14 specified on the boundary C of the domain is called 


a Dirichlet problem. If, instead of u, the normal 
derivative du/dn is specified on C, the problem 


is called a Neumann problem. 


DISCRETIZATION Discretization is the replacement of a derivative 
8: page 17 by a finite-difference formula (e.g., the FORWARD- 


DIFFERENCE FORMULA). 


DISPERSION The physical phenomenon of the velocity of a 

7: page 13 wave being dependent upon its frequency is called 
dispersion. 

DISPLACEMENT VECTOR The displacement vector Q0 is the correction 

ll: page 20 to gn at the (n+1)th iteration of an ITERATIVE 

$: page 78 SCHEME, i.e., 


ew - xr NMCE 


te 
DOMAIN (1) The domain of a given function f is/set of values 
for which f is defined. 
DOMAIN (2) A connected open subset of R” is called a domain. 
W: page 49 (A domain in one dimension is just an open interval.) 


DOMAIN OF DEPENDENCE The domain of dependence of the point (x,t) for a 
2: page 12 given HYPERBOLIC equation is the set of all points 
in the solution domain whose DOMAINS OF INFLUENCE 


include the point (x,t). 


DOMAIN OF INFLUENCE 
2 : page 12 


Wo: page 40 


EIGENFUNCTION 


EIGENVALUE 


pages 65,160 


Ux 


S : page 62 


ELLIPTIC 


EQUATION OF STATE 


l : page 15 


EVEN FUNCTION 


W : page 23 
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For a HYPERBOLIC equation, the two 
CHARACTERISTIC curves through the point (x,t) 
bound a domain which is called the domain of 
influence of the point (x,t). Values of the 
solution at x for a given time t can affect the 


solution only at points within this domain for 


which t > t. 
See EIGENVALUE. 


Given a linear problem consisting of the 

equation Lu + Apu = 0 (where L is a LINEAR 
OPERATOR) with homogeneous BOUNDARY CONDITIONS, 
it sometimes turns out that nontrivial solutions 
are possible only if A takes on prescribed values 
called eigenvalues. For a given problem, the 
solutions associated with each eigenvalue are 


called eigenfunctions (in the case of a function 


space) or eigenvectors. 
See CLASSIFICATION OF OPERATORS. 


The equation of state of a fluid is a relation 


between pressure and density. 
have the adiabatic equation pp Y = constant. 


A function f with domain R or [-a,a] or (-a,a) 


is even if, for each x in the domain, 


f(x) = f(-x). 


For an ideal gas, we 


EXPLICIT SCHEME 


`S : page |I 


FINITE-DIFFERENCE 
REPLACEMENT 


5 : page 11 


FINITE-DIFFERENCE 


SCHEME (or METHOD) 


5 : page 11 


FINITE FOURIER 
TRANSFORM 


W : page 129 


FINITE SINE 
TRANSFORM 


W : page 127 


FLUID EQUATIONS 


PDE Handbook 


A FINITE-DIFFERENCE SCHEME is explicit if each 
equation expresses one unknown PIVOTAL VALUE in 


terms of known pivotal values. 


A finite-difference replacement of a differential 


equation is a FINITE-DIFFERENCE SCHEME obtained by 


DISCRETIZATION of the derivatives in the equation. 


A finite-difference scheme is a formula relating the 


PIVOTAL VALUES of a function over a MESH. 


The finite Fourier transform of u(r,0), where 


-T < 0 < m, is the sequence la, G) b G2) given by 


T 
a(x) = las ng dg n= 


1 T E 
b, (1) 1 [ s¢e,o)sin n0 dð n= 


The finite sine transform of u(r,0), where 0 < 0 < m, 


is the sequence b (x)} given by 


2 


T 
b (r) = z u(r,0)sin n0 dO n= 1,2,.. 
n 0 


T 


See EQUATION OF STATE, MASS-CONSERVATION EQUATION and 


MOMENTUM EQUATION. 


FLUX 
3: page 21 


16: page 22 


FORWARD-DIFFERENCE 
OPERATOR 


5: page 9 


FOURIER-BESSEL SERIES 


FOURIER COEFFICIENT 


FOURIER SERIES 


W: page 72 
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The flux of a fluid is its rate of flow, i.e., 
the volume crossing a given surface per unit 
time. If the velocity of the fluid is given by 
the VECTOR FIELD v, the flux across the surface 


S is given by 


[rs ds 


where n is a unit vector normal to S. In the 
case of a fluid which flows down a pipe with 
speed w the flux through a cross-section D of 


the pipe is 


n w dA. 


The forward-difference operator ^n is defined by 


A, Y : x u(xth) - u(x) x € R. 


A Fourier-Bessel series is the FOURIER SERIES 
expansion of an arbitrary function with domain 


(0,1) as Y J CAS. in terms of BESSEL 


k=1 km 
Wem Am, _ 
FUNCTIONS satisfying VA ) =0. 


See FOURIER SERIES. 


The expansion of an arbitrary SQUARE INTEGRABLE 
function f as a series E ate. in terms of the 


ORTHOGONAL set of functions (9,7, with 
e = M 
“oe? 
k 9% 


is called a Fourier series with Fourier 


coefficients fe. 


GAMME FUNCTION 


14: page 7 


GAUSS-SEIDEL 
ITERATION MATRIX 
ll: page 23 


S: page 78 


GAUSS-SEIDEL 


METHOD 
li: pages 23,24 
S: page 78 


GLOBAL ERROR 


jun 


page 27 


loo 


page 7 


PDE Handbook 


The Gamma function is the function 


Tia — ee lay a> o 


For the ITERATIVE SCHEME 


xin) - ty 0D " yx +b, 


the matrix an lu is called the Gauss-Seidel 


iteration matrix. 


The Gauss-Seidel method for solving the problem 
(I-L-U)x = b, 

where L and U are respectively lower and upper 

triangular matrices, is given by the ITERATIVE 


SCHEME 


ux 


xr) = jy 0D * cobi 


m 
at a MESH point,[a FINITE- 


The global error, 
DIFFERENCE REPLACEMENT of a differential equation 
is the difference between the computed finite- 


difference solution and the true solution to the 


differential equation, at that point. 


GREEN'S FUNCTION FOR 


ORDINARY DIFFERENTIAL 


EQUATIONS 
W : page 122 
10: page 16 
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Green's function G(x,8) for the boundary value 
problem 


¿PY + q(x)u(x) = -£(x) xe (0,B) 
u(a) = u(B) = 0 
is given, for each E e (a,B), by the continuous 


solution of 


PoE] + acoc = 0 x4 E, 
lia z c| 790 

[egeo 7 S [eto 7 % 

dc _ dG 2.2.11 

dx x=E+0 dx x=E-0 JON 


Alternatively, we can write this system as 
d dG 
ax [POG] + ade = - 66-5, 


G(a,€) = G(B,E) = 0, 


where 6(x-§) is the DELTA FUNCTION in one 
dimension. The solution to the boundary value 


problem above may be expressed as 


B 
u(x) = Í G(x,8)£(5)d8. 


a 
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GREEN'S FUNCTION For any DOMAIN D with boundary C, Green's function 
FOR PARTIAL G(rip) for the boundary value problem 

DIFFERENTIAL Ltwl (x) z -F (x) reD 

EQUATIONS u(r) =0, rec 


Y : page 135 where L is a linear ELLIPTIC operator, is given, 


10: page 18 for each p e D, by the continuous solution of 
LEG] (£;p) = - §(r~p) re D 
G(x3p) = 0 re C. 
The solution to the problem above may be written 
in the form 
u(r) = | Gz39)F(g)dpo. 
D 
GRID See MESH. 
HARMONIC The nth harmonic of a function f with domain 
6 : page 14 [-L,L] is a sin mtx/L + b cos nxx/L, where a and b 
n n n n 
16: page 20 are the appropriate FOURIER COEFFICIENTS of f. 


HARMONIC FUNCTION A solution of Laplace's equation 


W : page 52 


Vu = 0 
is called a harmonic function. 
HEAT CONTENT Let u(x,y,z,t) be the temperature at the point 
3 : page 31 (x,y,z) in the domain D at time t > 0. The heat 


content of D is a measure of the total amount of 


heat in D and is piven by 


f 
[enne dxdydz 


JJ) 


D 
where p is the density and c the SPECIFIC HEAT 


of the material. 


HEAT EQUATION 


W: page 58 


HOLDER CONTINUOUS 


W: page 79 


HOMOGENEOUS EQUATION 


W: page 30 


HOOKE'S LAW 


W: page 6 


HYPERBOLIC 


IDEAL FLUID 


16: page 7 


IMPLICIT SCHEME 


S: page 18 


IMPROPER INTEGRAL 


6: page 14 


Glossary 


19 


The heat equation is another name for the 


DIFFUSION EQUATION. 


a 


The function f : I > R is Holder continuous 


at x € I if JM > 0, 0 > 0 such that Vy e I 


l£(y) - £(x)| s Mly-x!%. 


If L is any LINEAR operator then 


Lx = 0 


is a homogeneous equation. 


Hooke's Law states that force exerted by an 


elongated elastic medium is proportional to 


the extension per unit length; the constant of 


proportionality is called Young's Modulus. 


See CLASSIFICATION OF OPERATORS. 


An ideal fluid is one whose coefficient of 


VISCOSITY is zero. 


A FINITE-DIFFERENCE SCHEME is implicit if it 


must be solved for several unknown values 


simultaneously. 


A integral in which the integrand or the domain 


of integration is unbounded may be evaluated as 


an improper integral if it may be realized as 


the limit of a sequence of (proper) integrals. 


For example, if f(a*0) does not exist we may 


define 
m = lim | 
ja + jate 


provided this limit exists. 
ls . [b 
af = lim A f 
bm 


provided this limit exists 


Similarly we define 


20 


INDUCED INSTABILITY 


5 : page 6 


INFLUENCE FUNCTION 


W : page 119 


PDE Handbook 


À method of solution for a problem suffers from 
induced instability if the method magnifies LOCAL 
ERRORS so that the computed result differs 


significantly from the true result 


The influence function R(x,&) for the initial value 


problem 


pog] + g(x)u(x) = f(x)  x>0, 
u(a) = u' (a) = 0, 


is given, for each Ẹ > Q, by the solution of the 


initial value problem 


ipw] +q(x)R=0 x> &, 
R| ="05 
dE Bob. 
dx -E p(£) > 


or the equivalent problem 
dí, dR ur 
PO] + aGor = 56-5, 
R| ee EI 


The solution of the original problem is given by 


x 
u(x) = | R(x,5)£(6)d8. 
o 


INHERENT INSTABILITY 


5: page 5 


INITIAL-BOUNDARY 


VALUE PROBLEM 


INITIAL CONDITIONS 


l: page 9 


INITIAL VALUE PROBLEM 


ITERATIVE SCHEME 


ll: page 15 
JACOBIAN 
15: page 15 
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A mathematical problem suffers from inherent 
instability if it is not CONTINUOUS WITH 


RESPECT TO ITS DATA or is not WELL CONDITIONED. 


See CAUCHY PROBLEM. 


Initial conditions for a partial differential 
equation are conditions given at spatial points 
for time t - to: For example, a stretched 
string may have its motion started with a 
given shape, u(x,0) = f(x), and a given 
velocity distribution 356,0) - g(x). For 
ordinary differential equations, initial 


conditions are conditions specified at one point 


of the solution domain. 
See INFLUENCE FUNCTION. 


An iterative scheme for solving a problem is 
one in which the solution is obtained by a 


process of successive approximation. 


The Jacobian of the N functions Early 


: » r N 
each with domain variable (X 929 ++) eR 


and codomain R is the determinant of the matrix 


9f, af, | 
x, Ixy 


i, Oy 


22 


JACOBI METHOD 
ll: page 24 


S: page 78 


JUMP DISCONTINUITY 


6: page l4 


KINETIC ENERGY 


2: page 14 


KRONECKER DELTA 


LAPLACE'S EQUATION 


W: pages 43,49,52 


LAW OF CONSERVATION 
OF ENERGY 


2: page 13 


PDE Handbook 


The Jacobi method is the ITERATIVE SCHEME 
(n+) | 
x = 


for the solution of the equation 
(I-L-W)x = b, 


where L and U are respectively lower and upper 


triangular matrices. 


The function f has a jump discontinuity at x if 
the one-sided limits f(x+0) and f(x-0) exist 


but are not both equal to f(x). 


The kinetic energy of a particle is the energy 
it possesses by virtue of its motion, and is 


defined to be 


(mass) x (velocity)?. 


The Kronecker delta is given by 


fo i#j 


ij 1,j € Z. 


Laplace's equation is V*u = 0. 


The Law of Conservation of Energy states that the 


change in the total energy (= KINETIC ENERGY + 


POTENTIAL ENERGY) of a system is equal to the work 


done on the system by all the forces. In 
particular the energy of an isolated system 


remains constant. 


LAW OF CONSERVATION 
OF MASS 


l: page 12 


LINE INTEGRAL 


3: page 12 


LINEAR (OPERATOR) 


W: page 29 


LINEAR PARTIAL 


DIFFERENTIAL EQUATION 


W: page 30 
LOCAL ERROR 
5: page 27 


LOCAL TRUNCATION 
ERROR 


5: page 14 
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The Law of Conservation of Mass states that the 


total mass of a quantity of material remains 
constant in time, however it becomes distributed 
in space. 

f ds represents lim Xf, As 


C 
As-» 0 
over the CURVE C, where fp is the value of f 


The line integral | 


at a point P in the small element of C with 


length As. 


L is a linear transformation (operator) if it 
has the property 

Liou + By] = oL[u] + BLIV], 
where Q and B are any constants and u,y are any 


vectors (functions). 


Llu] = F 
where L is a LINEAR OPERATOR involving partial 


differentiation and F is a given function. 


The local error is the error made at a point 
when a FINITE-DIFFERENCE SCHEME is used once 
only. In the complete step-by-step process, 
the finite-difference scheme is used many times 
so that local errors accumulate to produce 


GLOBAL ERRORS. 


The local truncation error (at a point) in a 
FINITE-DIFFERENCE REPLACEMENT of a differential 
equation is the LOCAL ERROR introduced by 


DISCRETIZATION. 
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MASS-CONSERVATION The Mass-Conservation equation expresses 


EQUATION mathematically the physical LAW OF CONSERVATION OF 
l: page 13 MASS as applied to problems in fluid flow. It 
14: page 26 is given by 

30 A 

— + = 

3c div pu = 0, 


where p and u denote the density and velocity 
respectively, and t is the time coordinate. In 
one spatial dimension with coordinate x the 


equation reduces to 


9p du 90 _ 

at Pax * "ax ^O 
MEAN SQUARE The mean square deviation of g from f, with 
DEVIATION respect to the weight function p on [a,b], is 
W: page 71 


b (e Q2 
f (f-g) “p. 


MESH The set of points { (x; 8,3) = ((xptin,t +jk)} 


0 
5: page 7 where h and k are constants is called a mesh 
or grid. In a given mesh each point eut) is 


called a mesh point, and h and k are called the 


mesh spacings or mesh lengths in the x- and t- 


directions respectively. 


MESH RATIO _ A FINITE-DIFFERENCE REPLACEMENT of a HYPERBOLIC or 

5: pages 11,25 PARABOLIC differential equation depends on the 
MESH spacings h and k (in the x- and t-directions 
respectively) according to the value of the mesh 
ratio k?/h? (for hyperbolic equations) or k/h? 
(for parabolic equations). For typical examples, 


see Section 3. 


MOLECULAR DIAGRAM 


S: page 12 
5: page 11 
MOMENT 

16: page 13 


MOMENTUM EQUATION 
l: page 15 


14: page 27 


NEUMANN PROBLEM 


NEWTON 


NEWTONIAN FLUID 


16: page 7 
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A molecular diagram is a diagrammatic 
representation of a FINITE-DIFFERENCE SCHEME 
which indicates the coefficients of the PIVOTAL 


VALUES related by the scheme. 


The moment about the point P of a force, acting 
through the point Q and represented by the geometric 
vector F, is r x F where r is the geometric 


vector PQ. 


The Momentum equation expresses mathematically 
NEWTON'S SECOND LAW OF MOTION as applied to fluid 


flow. It is given by 


where p,u,p and F denote the density, velocity, 
pressure and BODY FORCE per unit mass respectively, 
and t is the time coordinate. In one spatial 
dimension with coordinate x the equation reduces 


to 


See DIRICHLET PROBLEM. 
A newton is the SI unit of force. 


A Newtonian fluid is one whose coefficient of 


VISCOSITY does not depend on the velocity gradient. 
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NEWTON'S SECOND LAN Newton's Second Law of Motion states that 


OF MOTION force = rate of change of linear momentum 

l: page 7 7 mass x acceleration (for a body of constant mass) 
NUMERICAL The numerical characteristics of a FINITE- 
CHARACTERISTICS DIFFERENCE SCHEME are the lines which bound the 

5: page 13 NUMERICAL DOMAIN OF INFLUENCE. 

NUMERICAL DOMAIN The numerical domain of dependence of a MESH 

OF DEPENDENCE point (x,t) for a given FINITE-DIFFERENCE SCHEME 

8: page 12 is the set of mesh points whose NUMERICAL DOMAINS 


OF INFLUENCE include (x,t). 


NUMERICAL DOMAIN The numerical domain of influence of a MESH point 
OF INFLUENCE (x,t) for a given FINITE-DIFFERENCE SCHEME consists 
3: page 12 of all those mesh points at which the PIVOTAL 


VALUES of the solution depend on the value of the 


Solution at (x,t). 


ODD FUNCTION A function f with domain R or [-a,a] or (-a,a) 
W: page 12 is odd if, for each x in the domain, 


f(-x) = -f(x). 


ORTHOGONAL Two elements, a and b of an inner product Space 
W: page 71 , are orthogonal if a'b = 0. For example, two real 


functions $ and Y are orthogonal with respect to the 
weight function p on (a,8) if 

B & 

[Povo = 0. 


ORTHOGONAL BASIS An orthogonal basis is a COMPLETE set of 


ORTHOGONAL functions. 


ORTHONORMAL FUNCTIONS 


6 : page 12 


OVER-RELAXATION 


FACTOR 


PARABOLIC 


PERIODIC FUNCTION 


PIECEWISE CONTINUOUS 


6 : page 15 


PIECEWISE CONTINUOUSLY 


DIFFERENTIABLE CURVE 


PIVOTAL VALUE 


5 : page 7 


POINTWISE CONVERGENCE 
(OF A SERIES) 


W : page 70 
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A set of functions [e which satisfies 


8 
ain), 090 dx = Spn? 


where [aM is the KRONECKER DELTA, is said to be 


orthonormal with respect to the weight function 


P on the interval (a,8). 


See SUCCESSIVE OVER-RELAXATION METHOD. 


See CLASSIFICATION OF EQUATIONS. 


A function f is periodic (with period a)if, for all 


x € R, f(x) = £(xta). 


A function is piecevise continuous on [a,bl if it is 
continuous at each point of [a,b] except for 


a finite number of JUMP DISCONTINUITIES. 


See CURVE. 


The values of a function u at MESH points are 


called pivotal values of u. 


Given an infinite sequence of functions 
9,595, with domain I, we let 
N 
= . I 
Sy (x) B=] Enn 09 If for each x € 


lim Sy 6) = f(x), we say that the series 
N yo 


Pe 00 converges to f(x) pointwise in I. 
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POISE A poise is the cgs unit of VISCOSITY. 


16: page 6 


r 3 i E 
POISSON'S EQUATION Poisson's equation 1s 


ia 

W : pages 49,50,129 Vu = -F 
3 : page 5 
POTENTIAL ENERGY The potential energy of a system is the energy 
2 : page 14 stored within the system. 
PROPERLY POSED A system consisting of a partial differential 
PROBLEM equation together with (INITIAL and) BOUNDARY 
W : page 6 CONDITIONS is properly posed if its solution 
3 : page 10 (a) exists, (b) is unique and (c) is CONTINUOUS 

WITH RESPECT TO THE DATA. 
PROPERTY A À square matrix is said to have Property A if, 
ll: page 29 by transposing pairs of rows and pairs of 


corresponding columns, it can be transformed to 


the form of the matrix 


where bi and D, are diagonal matrices. 


RAYLEIGH QUOTIENT Consider the differential equation 


W : pages 163,165 (pu')' - qu * Apu = O in. (0,8), 


where q and p are continuous and p is CONTINUOUSLY 
DIFFERENTIABLE in (0,8), along with u(a) = u(B) = 0. 
Let $ vanish at Q and B and be twice continuously 
differentiable in (a,B). The Rayleigh quotient of 


$ is given by the ratio 


RECURRENCE RELATION 


ROUNDING ERROR 


RESIDUAL VECTOR 
11: page 16 


S: pages 31,79 


SCALAR FIELD 


3: page 6 


SELF-ADJOINT FORM 


W: page 117 
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Be ards 42 
[Bove *qQ^) 
Boa 
[epe 
Characterizations of the EIGENVALUES of the 
problem above as minima of the Rayleigh 


quotient are called minimum principles for the 


eigenvalues. 


À recurrence relation is an equation relating 


Successive members of a sequence. 


A rounding error is introduced into the solution 
of a problem by working to a definite number of 


significant figures. 


Given a system of linear equations 
Ax = b, 
then, if x is an approximation to x, the vector 


E-b-Ax 


is called the residual vector. 


Let D denote a subset of 83; then a function 
$:D^R $:Dx RJ, +R 
or 
$ : (x,y,z) u $ : (x,y,z,t) u 


is called a scalar field. (In the second case 


we say that $ is time-dependent.) 


An ordinary differential equation written in the 
form 

(pu')' + qu - f 
on the interval (a,b) is said to be in self- 
adjoint form. We require that p be CONTINUOUSLY 
DIFFERENTIABLE and positive and that q and f be 


continuous on [a,b] 
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SINGULAR POINT 


10: page 7 


SOR 
SPARSE 


ll: page 15 


SPECIFIC HEAT 


3 : page 31 


SPECTRAL RADIUS 
8 : page 18 


S : page 79 


SQUARE INTEGRABLE 


6 : page 13 
STABLE 
8 : page 17 


A singular point of the ordinary differential 


equation 
a) + NES + c(x)u = 0 
dx dx 
is a point x where either a(x) = 0 or b(x) or c(x) 
is infinite. 
See SUCCESSIVE OVER-RELAXATION METHOD. 


A matrix is sparse if it has a large proportion of 


zero elements. 


The specific heat of a material is the quantity of 
heat required to raise the temperature of unit mass 


of the material by one unit. 


The spectral radius of the matrix M is 


p™) = maxla; | 
i 
where A, G=1,2,... ,n)are the EIGENVALUES of M. 


The function f is square integrable on (es) with 


respect to the weight function p if the integral 


B 
£?p 
0 


exists. 


A FINITE-DIFFERENCE SCHEME is stable if LOCAL ERRORS 
do not accumulate exponentially as the step-by-step 
solution progresses, i.e.,if it does not suffer 


from INDUCED INSTABILITY. 


STANDARD FORM The standard form (or canonical form) for a 


W: pages 42,43 linear partial differential operator which is 


.232u du du 
(a) HYPERBOLIC is Eon * FGp pu Em = 05 
22 9?u du du 
(b) PARABOLIC is »w * FGp pu Em = 05 
(c) ELLIPTIC is uy, Lam reu, om -0 
dE an de sn.” : 
STANDING WAVE A standing wave is a wave which does not move 
1 : page 35 laterally, but merely changes its amplitude 
2 : page 9 with time. A standing wave is easily recognized 


by observing points of zero displacement which 


remain stationary for all time,e.g., Asin wt cos kx. 


SUBSIDIARY CONDITIONS Subsidiary conditions are INITIAL and/or 


BOUNDARY CONDITIONS. 


SUCCESSIVE OVER- The successive over-relaxation method for 
RELAXATION METHOD solving 

(SOR) (I-L-U)x = b, 

ll: page 26 where L and U are respectively lower and upper 


triangular matrices, is given by the ITERATIVE 
SCHEME 


(n) 


(uy = {(1-w)I + wh"? + wb. 


The parameter œ is called the over-relaxation 
factor, and determines the ASYMPTOTIC RATE 
OF CONVERGENCE of the method. The optimum 
over-relaxation factor is that value which 


maximizes the rate of convergence of the method. 
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SURFACE INTEGRAL The surface integral E dA of the function f 


3: pages 11,27 over the surface S represents lim Zf AA, where 
AAO 
fp is the value of f at a point P in the small 


element of S with area AA. 


TAYLOR APPROXIMATION A Taylor approximation of a function is obtained 


l: page 13 by neglecting the remainder term in TAYLOR'S 
2: page 21 THEOREM (see Section 3). The first-order Taylor 
14: page 26 approximations for functions of one,two and three 


variables are respectively: 
u(xth) = u(x) + hu! (x); 
du du 
ð 


u(x+th,t+h) = E + hx * ke (x, t)! 


u(x+Ax,y+Ay,z+Az) = [u + Ar*grad Ud 353) 
> , 
where Ar = (Ax,Ay,Az). 


THERMAL CONDUCTIVITY The thermal conductivity of a material is the 


W: page 58 rate at which heat flows across a temperature 


gradient of unit magnitude. 


THREE-LEVEL SCHEME A three-level FINITE-DIFFERENCE SCHEME involves 
15: page 6 terms along the (j-1)th, jth and (j+1)th time 
levels. If it is linear, it may be written in 


matrix form as 


TRAVELLING WAVE 


l: page 20 


TRIDIAGONAL MATRIX 


5: page 18 


UNIFORMLY BOUNDED 


6: page 21 


UNIFORM CONVERGENCE 


2: page 10 
6: page 9 
W: page 70 
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A travelling wave is a wave which moves in a 

given direction without changing its shape. 

For example f(x*ct) represents’ a wave 

travelling in the negative x-direction. The 
Substitution X - x*ct gives f(X) which may 
represent a wave of fixed shape. Since the origin 
of the X-coordinate is at x - -ct, the wave 
travels with velocity c in the negative x- 
direction. Similarly, f(x-ct) represents a 
travelling wave of velocity c in the positive x- 


direction. 


A tridiagonal matrix is one whose nonzero elements 
appear on and adjacent to its main diagonal, i.e., 


a BAND MATRIX of BAND WIDTH 3. 


The infinite sequence of functions (6) is 
uniformly bounded on the interval I if Hc e R, 
independent of n e z* and x e I, such that 


ne Z' and x c T4 lb 09] * c. 


The infinite sequence of functions (£) 
converges uniformly to f on the interval I if 
it is convergent to f in the uniform norm, 


È 
i.e., if for each € > 0 dw, € Z such that 


max |f(x) - f£ (x)| <e 
N 
xel 


whenever N > Ne. An equivalent formulation is: 
for each € > 0,3 No e z* independent of x such 
that |£G0 - Ey 09 | < € for all x e I whenever 


N2N. 
E 
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VECTOR FIELD 


3: page 6 


VELOCITY POTENTIAL 


14: page 15 
VISCOSITY 
3: page 20 


16: page 6 


PDE Handbook 


E 
The infinite series of functions DECLA 


converges uniformly on the interval I if the 


sequence {s } of partial sums 


is uniformly convergent. 


A vector field in two (or three) dimensions is 

a function with domain D in R? (or R?) and codomain 
the space c? (or G3) of geometric vectors. A time- 
dependent vector field is defined similarly with 


de 
domain D x Ro: 


If the velocity y of a fluid is related to a 
SCALAR FIELD 9 by v = -gradó, 9 is called the 
velocity potential of the flow, which is then said 


to be irrotational. 


If a relative motion occurs in a fluid, a 
measurable resistance is experienced, and the 
fluid is said to exhibit viscosity or internal 
friction. For a fluid moving in parallel layers, 
the frictional force F on a plane surface of area 
A parallel to the fluid flow is given to a good 


approximation by 


where Ss is the velocity gradient perpendicular 
to the direction of motion and u is called the 


coefficient of viscosity. If p is the density of 


the fluid, then its kinematic viscosity is p/p. 
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WAVE EQUATION The wave equation is 
W : pages 9,36,48,66,152 3?u 292 
—— = c Vu 

2 
l : page 8 at 
WAVE LENGTH See ANGULAR FREQUENCY . 
WAVE NUMBER See ANGULAR FREQUENCY. 
WELL CONDITIONED The linear equation 
ll: page 15 Au = b 


is well conditioned if small changes in the 
elements of A and b do not cause large changes 


in the solution vector u. 


YOUNG'S MODULUS See HOOKE'S LAW. 
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2 VECTORS 


2 
A vector space | over a field F of sealars is a set of elements with an internal binary 
operation &£:E » 4 — E and an external binary operation of. multiplication 
F x V —À9 V satisfying the following axioms for all a. b. e« Fand mone F: 
| (a+b) +e =a + (b+ ce) 

2 a+b=b+a 
3 30eV such thu 0 a za. Vae V 
4 3-aeV such thata + (7a) — 0 


5 mía b) = ma + mb 


6 (m na = ma + m 
7 (mnja = mina) 
$ la=a 
When F = R (the reals) we talk about a real rector space: when F = C tthe comples 
numbers) we have a complex vector space. 
An important notion in a vector space is that of linear independence. A set !v,..... Yu 
of vectors in V is linearly independent if 

AV T Xa bosco xw, =0 (z,€ F) 


holds only for z, = 72 27: = z, - O Suppose now that fy,... Yas isa maximal 


ans 
linearly independent set in V. i.c, if another vector in V is included in the set we obtain a 
set which is linearly dependent. Then [v,..... Ym, i5 called a basis for V. and every 
other basis for V has precisely m vectors. We say that V has dimension ni. It can also 
be shown that an arbitrary element x ¢ V can be expressed uniquely in the form 


NENN tot X, (x, € F). 


— 
Geometric vectors are equivalence classes of parallel arrows of equal length. If AB is 
an arrow, we write AB for the geometric vector to which it belongs. The sum of two 
geometric vectors is obtained by choosing suitable arrows which can be added accord- 
ing to the triangle law of addition: 


OA + AB = OB. 
So o > 


QNT. LS ^ 
If 4 is a real number then we define the product ¿AB to he the geometric vector 
whose arrows are parallel to those of AB but with lengths multiplied by 7 


[e] 


With respect to these definitions of addition and scalar multiplication the set G of 
geometric vectors in space forms a real vector space. We often consider the set Got 
geometric vectors in the plane, which is a subspace of this vector space. The null 
vector 0 is the geometric vector whose arrows are all of zero length. 


Let Oxyz be a frame of rectangular Cartesian coordinate axes and let i, j. k be 
geometric vectors whose arrows are of unit length and parallel to the coordinate axes 
Ox. Oy. Oz respectively. Then every geometric vector in the plane can be written 
uniquely in the form 


a = ai + aaj. 
and every geometric vector in space can be written uniquely in the form 
a =i t hj + ayk 


Thus ʻi. j! forms a basis for the space G? of geometric vectors in the plane, and ji. j. kj 
farms a basis for Gy 4. 05. 04 are, respectively. the v-. z-coordinates of a. 
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If I is an interval on the real line or a domain in R? or R?, then the 
set of all functions I + R is a real vector space with addition and scalar 


multiplication of functions given by 


(f+g) : xı f(x) + g(x) x €t, 


Vie R (Af) : x PALE (x)] xe I. 


Replacing R by C in the above definition yields a complex vector space 


of functions. 


The space of all VECTOR FIELDS with domain D forms a real vector space 


with addition and scalar multiplication of fields given by 


(uty) : ch ulr) + v(x) x € D, 


NACER Ay : re AyG)l . rcD. 


An important property of à geometric vector is its lengrh or magnitude. To define this 
more generally, we introduce the inner product. 
A mapping *: V « -=> R. where Eis a real vector space. is called a (reall inner 
product if for alla. b.e Land ^g R: 
| aca sera 0: 
2 ara? 0: 
5 aeb bea: 
dl o adeb os geb ra c paee 
We mias specily an inner product on G^ or G* by 
OA-OB |. OAXOBxvos 0, 
— — . Ly . . 
where OA. OB denotes the length of the arrow OA, OB respectively and 0 is the angle 


Ls: d —Ó 
between OA and OB. The geometric vectors i, j. k are orthogonal, Le. 


iej- jeko kei -0. 


Un fact any two geometric vectors which are represented by perpendicular arrows are 
T 3 E . M A H a k 
orthogonal.) Tf OA= eX * Say + ask cM OB = Et ba) th h, thom 


OA ORB ahy toba 9 au. 
in terms of coordinates 


Given the vector space of SQUARE INTEGRABLE functions with domain (0,8) 


and codomain R, the map 


8 
(f,g)h— | fg 
a 


specifies an inner product. 
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A mapping * : V x V + C, where V is a complex vector space, is called a 
complex inner product if for all a,b,c e V and A,u e C: 
| ara=0é5a=0; 


2 ata 


E 


3 2b=bra 


4  g'Ob*uc) = larb + Marc. 


Schwarz's Inequality states that, given a real inner product space V, 


1 1 
xy S EDEN” 
for all x,y € V. For the space of square integrable functions with the 


inner product given previously, this becomes 
fg 8,245 fB 233 
i < 2 212 
Beg s Ji 


A mapping || ||: V > R, where V is a real vector space, is called a 


norm if for all a,b e V and A e R: 


1 all = 0 >a = 0; 

2 a|| 20; 

3 [Dell = I3 lla ll 

4 arbli < Ilall +11 bl | (Triangle Inequality). 


In an inner product space the specification 


llall= (ara) 


yields a norm satisfying the axioms above. The norm of a geometric 
vector, given in this way, is just the length of an arrow representing 
it; by convention, if a represents a geometric vector (or a vector field) 


then its length is denoted by la] or just a. 
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In G* pbut not in G% se define the vector product x 360. G* - + GS by 


OX x OB ON OB sin ten, 
=> => 


> -> 


where fis theanglet 5 0 «. zi between OA and OB, and n is the unit rector (Le. vector 
Whose norm is 1) orthogonal to OA and OB, such that OA. OB and n form a right- 


banded triad. 


The vector product has the following properties: 


1 
" 
A 
3 


n 


ixjekjxk -ikxi- jixisjxjekxk -0: 


axb= -bxa: 


axib-c)j-axb-axce: 


a x b = (usb, = aub = tahy 


"bj ta hs = aha: 


(ax b) xe 4a x (bx c). in general. 


E 
ES 
ar 
x 
[e] 
2, 
" 
c 
e 
PA 
io 
x 
in 
2 
" 
A 
A 
v 
x 
e 
Z 
[ 


8p 9$ Si 
ay Dc 
a, b3 93 


A point in space with Cartesian coordinates (x,y,z) may also be 


specified by its spherical polar coordinates (r,0,0) or by its 


cylindrical polar coordinates (0,0,2). These are related to each other 


as follows. 


r? = x2 y? + 2? 
x = r sin O0 cos bd 
y= r sin 6 sin > 
z= r cos 0 


2 2 


oy. 


p = x 
b 


p cos 


p sin $ 


In evaluating a VOLUME INTEGRAL the volume element dV is replaced by 


dx dy dz 


r?sin 0 dr de dd 


p dp do dz 


in Cartesian coordinates 
in spherical polar coordinates 


in cylindrical polar coordinates 
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The slope of a SCALAR FIELD Y is given by the VECTOR FIELD 


A, ay 
Brad Y = gaiti eE 
LL 1 Y, 


where EN denotes the vector field of unit norm in the direction of 


increasing a. 


If v is a VECTOR FIELD given by 


y^vi*vi*vk 
^ x^ z 


=y * Vg&g * VEZ 


8 $76 


=ve +ve,+vk, 
pap $76 z^ 


e 
ETE 


we may define the SCALAR FIELD 


| 
| 


div v 


The LAPLACIAN operator is given by 


y?y 


div grad Y 


2 2 2 
ay mu. an 
ax? ay? 3z? 
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The following rules are useful: 
grad(¥) = Y grad 9 + È grad Y 
div(Yv) = Y div y + v'grad Y 


div(Y grad 9) = Y V?%9 + grad V'grad 6 


Divergence Theorem (3: pages 14,28) 
Let D be a CONVEX DOMAIN in R? (or R?) bounded by the closed CURVE (or 
closed surface) C. If the coordinates of the VECTOR FIELD v are 


continuous on DuC and their first partial derivatives are continuous in 


D then 
f i [ 2 
] iv y dà = Ó ven ds (D c RÊ), 
D J 
f div y dV = | vin d$ (D c R°), 
Jp C 


where n is the unit outward normal to C. 


Green's Theorem (W: page 53) 
Let D be a CONVEX DOMAIN in R? (or R?) bounded by the closed CURVE (or 
closed surface) C. If the SCALAR FIELD u has continuous first partial 


derivatives on DuC and continuous second partial derivatives in D then 


u V?u dA = j u X ds - | lgrad u|?dA (D c R2), 
lp e D 
f u Vu dV = [ u ou ds - | | grad u|?av (Dc R’); 
Jp 4g * D 


where du/dn is the outward normal derivative of u on C. 
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3 USEFUL FORMULAS AND THEOREMS 


Trigonometric Functions 


sing - 8) = cosé 


cos G - 8) sin8 


sin?8 + cos?0 = | 
sin20 = 2sin@cos® 
cos20 = cos?0 - sin?0 
= 2cos?0 - 1 = 1 - 2sin?8 
sin30 = 3sin9 - 4sin?0 
cos30 - 4cos?0 - 3cos0 


tan20 = 2 tanð — 


1 - tan?8 


1 + tan?@ = sec?8 


1 + cot?@ = cosec?8 


sin(atB) = sinocosB + cosasinf 

cos(atB) = cosacosfí + sinasinf 
tana + tan 

tan(atB) = tano + tang 


1 + tangtang 


sina + sinB = 2sin }(a+tB)cos}(a-8) 


sina - sing = 2cos l(o*f)sin 3(a-B) 


cosa + cos8 = 2cos l(a*8)cos }(a-f) 


Li 


cosa - cos = -2sin l(os8)sin l(a-B) 
cotO + tan@ = 2cosec20 
cot8 - tan8 = 2cot20 


¿e cos® + isin (Euler's formula ) 
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Hyperbolic Functions 


: x — : . . 
sinh x = l(e" - e ^) = -i sin ix 
x “Xx . 
cosh x = l(e^ + e ^) = cos ix 
sinh x 
tanh x = 
cosh x 


cosh?x - sinh?x = | 


Trigonometric Fourier Series 


The trigonometric Fourier series of the function f with domain [-L,L] 


is given by 


oo 
nix + nx 
f(x) v 22, + a (s, ces mp b sin em) 5 
where 
1 [L nTx 
arr [E #0008 dx 
A Jb + NTX 
ba —T fi £cosin L dx 


Taylor's Theorem (5: page 8, 8: page 10) 
If u is a function of one variable which is n+] times continuously 
differentiable on the interval [x,x+h], then there is a 0 e (0,1) such 


that 
n+] 


DT ¿e (x*0h). 


n 
"u(Gxh) = ux) + hu'(x) +... + Bu qo + 


m à dass +1), 
(for h < 0, [x,x*h! is replaced by [x*h,x].) If, in addition, yA ) is 


bounded on [x,xth], i.e., 


pe) Gy) eB x Daan] 


j 3 n+l 
for some B, then the remainder term is of order h A 


antl 


uer uU P Geom) = omt!) as nago. 
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Intermediate Value Theorem (8: page 10) 


If a function f is differentiable on the interval [a,b], then for any 
real number y between f'(a) and £'(b), there is at least one point 


xy € [a,b] such that £' (xp) = y. 


0 


Principle of Superposition (W: page 33) 
The solution of the linear equation 

L[u] = F 
subject to the linear conditions 

L;[u] = f; i = 1,2,...,n 


is given by 


ae a ae suce cde di. 
where uo satisfies the problem 
L[ugl =F 
L;Cug] = 0 x9 
and (for k = 1,...,n) "ue satisfies the problem 
Llu, J =0 


Liu = fiig 


(6 is the KRONECKER DELTA). 


ik 


Extremum Principles for Flux (3: pages 22,23) 


Let the SCALAR FIELD w satisfy V%w = -k in a domain D, with w = 0 
on C (the closed boundary of D). Then if w* is any nonzero SCALAR FIELD 
which is differentiable on D and continuous on DUC such that w* = 0 on 


C, and y* is any VECTOR FIELD such that div y* = -k in D, then 


el Jun a) 


s | w dA s — 
| Igrad w*]? dA D 
D 


and Theorems 


Maximum Principle for Poisson's Equation (W: 


If V%u = -F in D, and F < 0, then 


max u(r) < max u(r) 
reD reC 


where C is the boundary of D. 


Finite-Difference Formulas 


ul (x) = 5 CuGx) - uGcb)J 
1 


u' (x) = 2k Cu(xth) - u(x-h)] 
ul (x) = 2 CuGeth) - u(x) J 

2 2 
Bu p au, Ct, Oe 


9x?  83y? (Ax)? (Ay)? 


pages 55,56) 


45 


backward-difference formula 


central-difference formula 


forward-difference formula 


five-point formula 


(6. is the CENTRAL-DIFFERENCE OPERATOR) 


o 
N 
E 
^ 
* 
2 
l 


200 E 
h^u"(x) + 12 


The DIFFUSION OPERATOR 


au _ Pu 
dt ax? 


may be replaced by: 
the Crank-Nicolson scheme 
a = 
ga atra 


the Du Fort and Frankel scheme 


Yi j+ ^ Pagar 7 


u(xth) - 2u(x) + u(x-h) 


nia) ew) + O(h). 
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or the explicit scheme 


iz) = ?u, wa 
uijt “i,j o a 


Here r = k/h? is the MESH RATIO, and u, i denotes u(ih,jk). 
Courant-Friedrichs-Lewy (C.F.L.) Condition (8: page 14) 

A necessary and sufficient condition for the convergence of a FINITE- 
DIFFERENCE SCHEME for a HYPERBOLIC equation is that the NUMERICAL 
DOMAIN OF DEPENDENCE of a point must include the DOMAIN OF DEPENDENCE 


of the same point. 


Lax's Theorem (8: page 16) 
Given a PROPERLY POSED, LINEAR, INITIAL VALUE PROBLEM and a linear 
FINITE-DIFFERENCE REPLACEMENT which is COMPATIBLE with it, the finite- 


difference scheme is CONVERGENT if it is STABLE. 


A linear FINITE-DIFFERENCE REPLACEMENT of an INITIAL VALUE PROBLEM is 

STABLE if 

(i) the SPECTRAL RADIUS of the matrix corresponding to the 
finite-difference scheme is not greater than 1 (Matrix Method) 

or 

(ii) the solutions to the scheme which have the form ei bPh;a 


satisfy || < 1 (von Neumann's method). 
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Newton's Method 
The nonlinear system of equations 
E(x) =0, 
where f(x) denotes (£, 69.2, G0... £ G0), may be solved by the linear 


ITERATIVE SCHEME 


facta -0t gang 


where J (n) denotes the matrix whose determinant is the JACOBIAN of the 
x 

functions te) at x. 

Chain Rule 


If UX] Xo 00. Xy are variables which depend on t, and 


u= PX) 
then 
du. P ag 9X 
dt. , Ox. dt 
i=1 i 


Fourier Series 
Let to? be an ORTHOGONAL set and f SQUARE INTEGRABLE on [a,8] with 
respect to the weight function p, and let the Fourier series for f be 


9 
given by 


eo 
Then X c? [Bozo converges and satisfies Bessel's Inequality 
n 


n 
n=1 
» 2 8 2 8 2 
E € o*p < f^p 
n=1 "Jo P” a 


Equality holds if and only if the Fourier series for f CONVERGES IN THE 


MEAN TO f, in which case we haveParseval's Equation 
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If (e) is COMPLETE, and 
co 
ka Lckó, 
A apse 


then 


In the following, let L be the operator defined by 
L : ur (pu')' - qu, 
and let uk denote the kth EIGENFUNCTION of 


Lu + Apu = 0 in (a,B) 


u(a) = u(B) = 0. 


Lagrange's Identity 


[ce - vLu) = [p(uv' - vu!)1 


over any interval in which u and v are twice CONTINUOUSLY DIFFERENTIABLE. 


Theorem A 
The set tu is COMPLETE for the space of functions SQUARE INTEGRABLE 


on (a,8) with weight function p. 


Theorem B 


If f is continuous on [a,81, f(a) = £(B) = O and 


Bing? 2 
[fos + q£%) 


exists, the Fourier series of f in terms of lu) is UNIFORMLY 


CONVERGENT to f on [a,81. 
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Oscillation Theorem 


The eigenfunction u, has precisely k-1 zeros in the open interval (a,8), 


k 


and the zeros of successive eigenfunctions interlace. 


Tests for Uniform Convergence 


CAUCHY'S tls! 
A necessary and suficient condition for a sequence of functions u, with domain / 
to be uniformly convergent is that given any & > O there exists an integer N, inde- 
pendent of xe / such that 

lux) — ux x E 


for all i.i > N, and all x & Z. 


WHIERSTRASS’S M-1151 


‘This isa sulligient though not necessary condition, If each term of the series 


J= Ea) vel 


ned 
is positive and 
ma < My WrelneZ*. 


where each M, is independent of x. and if 


g 


Y M, 


wet 


is convergent, then the given series for f is uniformly convergent. 


Bessel Functions 
n is a solution of BESSEL'S EQUATION OF ORDER m, and is bounded at the 


origin. 
m+2k o 


; CO'Qo 


XS = Poo kino] 


L fI (e)] = e" J 
m m- 


) ) 


42 ICONE 
)'"^ where A LM ) 0 


(m) 
x3, CAU tax = ATT, UAE" 


Ji) = MUI, Ct) - J, E] 


50 


Function 


Gra)” 


2xta 


x^*axtb 


2xta 


(x? 


1 
Y (a?-x?) 


ERE NONE 
Y (a? ex?) 
1 
Y (x?-a?) 

sin x 
cos X 
tan x 
cot x 
sec x 
cosec x 
sectx 

a 
cosec x 


ax 
e 


+axtb)™ 
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TABLE OF INDEFINITE INTEGRALS 


Indefinite Integral with respect to x 


+1 

Ga)” - 

nt] (n # -1) 
£n|x-a| 
= tan! = 

a 

1 atx 

2lal is eas 


&n |x? +ax-4b | 


(x2+ax+b) 2 
I-n 


(n#1) 


=] 
+ cosh Al (sign that of x) 


sinh x 
cosh x 
tanh x 
coth x 
sech x 
cosech x 
sech? x 
cosech? x 


sin rx cos nx (r?tn?) 
sin rx sin nx (r? n?) 

24,2 
cos rx cos nx (r7#n*) 


ax 
e" cos(bx + c) 


a : 
e Xsin(bx * c) 
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cosh x 
sinh x 

£n (cosh x) 
2n|sinh x| 
2 tan`! (e*) 


2n | tanh ix| 


tanh x 
-coth x 
. cos(r-n)x _ cos(r*n)x 
2(r-n) 2 (rtn) 
sin(r-n)x _ sin(r+n)x 
2(r-n) 2(r*n) 
sin(r-n)x sin(r+n)x 
2 (r-n) 2(r+n) 
¿ex 
[acos (bx+c) + bsin(bx+c)] 
a?+p? 
¿ax 


Casin(bxtc) - bcos (bx+c)] 
a? +b? 
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In the following m and n are integers. 


le sin mx sinnx dx = 0 


T 
nxdx = 


nx 


nx 


2 

cos nx dx =(T 
J9 

cos nx dx = 


cos nx = 0 


dx = 


#0 


